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ABSTRACT 

We detect angular galaxy-QSO cross-correlations between the APM Galaxy 
Catalogue and a preliminary release (consisting of roughly half of the anticipated 
final catalogue) of the Hamburg-ESO Catalogue of Bright QSOs as a function 
of source QSO redshift using multiple cross-correlation estimators. Each of the 
estimators yield very similar results, implying that the APM catalogue and the 
Hamburg-ESO survey are both fair samples of the respective true galaxy and 
QSO populations. Though the signal matches the expectations of gravitational 
lensing qualitatively, the strength of the measured cross-correlation signal is sig- 
nificantly greater than the CDM models of lensing by large scale structure would 
suggest. This same disagreement between models and observation has been found 
in several earlier studies. We estimate our confidence in the correlation detections 
versus redshift by generating 1000 random realizations of the Hamburg-ESO QSO 
survey: We detect physical associations between galaxies and low-redshift QSOs 
at 99% confidence and detect lensing associations at roughly 95% confidence 
for QSOs with redshifts between 0.6 and 1. Control cross-correlations between 
Galactic stars and QSOs show no signal. Finally, the overdensities (underdensi- 
ties) of galaxies near QSO positions relative to those lying roughly 135 — 150' away 
are uncorrelated with differences in Galactic extinction between the two regions, 
implying that Galactic dust is not significantly affecting the QSO sample. 

Subject headings: gravitational lensing: galaxies: catalogs: large scale structure 
of the Universe 
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1. Introduction 

Turner (1980), noting the 'apparent extreme evolution' of the QSO luminosity function, 
first indicated along with Canizares (1981) that magnification bias due to gravitational 
lensing could create statistical associations between foreground galaxies and high-redshift 
QSOs. A number of studies have since examined the phenomenon. 

Magnification bias arises when foreground gravitational lenses (i.e. galaxy clusters or 
large scale structure) magnify the source planes of background objects (QSOs). This source 
plane magnification results in two effects: 1) Individual sources span a larger area on the 
sky than they would if a lens were not present. Because gravitational lensing conserves flux, 
this results in a net brightening of the source, which can lead to the inclusion of otherwise 
undetectable objects in a flux-hmited sample. 2) The number density of the sources on the 
sky is decreased because of the same stretching of the source plane by the lens. 

Assuming a flux-limited sample of QSOs, depending on the slope of the magnitude num- 
ber counts of a sample of selected sources, am, the combination of the two effects produces a 
statistical association, ujoqid), or alternatively an overdensity (or underdensity) , q{0), which 
is dependent on the magnification, fJ,{9): 

q{e) - 1 = ugq{B) = iJiief - 1. (1) 

If am > 0.4, positive correlations will be observed, whereas negative correlations will occur if 
am < 0.4. Should am equal the critical value of 0.4, no correlations will be detected regardless 
of lens magnification. Finally, the amount of magnification depends on the projected surface 
mass density of the lens and on the relative distances of sources and lenses from the observer. 

There are, however, complications. For instance, if dust were located inside galaxy 
clusters, then QSOs lying on lines of sight near foreground clusters would tend to be more 
heavily obscured than those not lying behind clusters. This would result in a net decrease 
in a measured galaxy-QSO cross-correlation. A wide variety of techniques have been used 
to search for dust in galaxy clusters spanning their range in cluster richness. For a more 
complete description of past work in this area, please refer to the introduction section of 
NoUenberg et al. (2003). Early studies performed on the basis of visual inspection (Zwicky 
1962; Karachentsev & Lipovetskii 1969; de Vaucouleur, de Vaucouleur & Corwin 1972) 
tended to report rather large amounts of extinction (up to Ay ~ 0.4 magnitudes) in clusters 
such as Coma and Virgo. Smaller (though occasionally mutually inconsistent) estimates 
have been reported by later studies (Wise et al. 1993; Hu 1992; Bahcall 1999; Arnaud 
& Mushotzky 1998; Voit & Donahue 1995; Annis & Jewitt 1993; Romani & Maoz 1992; 
Veron-Cetty & Veron 1989; Maoz 1995; Dwek, Rephaeh & Mather 1990; Ferguson 1993; 
Biviano et al. 1990; NoUenberg et al. 2003), including upper limits of E{B -V) ^ 07^02 
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and Ar — CP025 on reddening and extinction by large scale distributions of dust within 
medium- richness APM clusters, respectively (Nollenberg et al. 2003). Low levels of dust 
in galaxy clusters such as this imply that cluster dust does not significantly contribute to 
galaxy-QSO associations, though the matter has not yet been entirely resolved because the 
possibility exists that dust may reside in small clumps which could have escaped detection 
using the Nollenberg et al. (2003) technique. 

Galactic dust is also a thorny issue because its presence could result in positive cor- 
relations between galaxies and QSOs or between galaxy clusters and QSOs. Due to the 
fact that Galactic dust extinguishes light from galaxies and QSOs alike, the observer would 
tend to see QSOs and galaxies that lie in the same unobscured parts of the sky, yielding a 
positive statistical association. Its distribution is also quite filamentary and this structure 
continues to scales smaller than are sampled by either the ~ 6'1 pixel size of the Schlegel, 
Finkbeiner & Davis (1998) extinction map or the lower resolution Burstein & Heiles (1982) 
map. High- contrast filaments of dust could therefore obscure individual sources, and while 
statistical association measurements generally sample scales much larger than this. Galactic 
dust could in principle leave its imprint by introducing power on a wide variety of scales that 
are dependent on the topology of the dust distribution itself. 

This paper is organized as follows: After a discussion of prior results in Section 2, we 
will describe the galaxy and QSO data that we use in Section 3. Then we will discuss 
our methods for estimating galaxy-QSO cross-corrclationsin Section 4. Following this, we 
will discuss the results of our analysis, including measurements of angular cross-correlations, 
integrated cross-correlations vs. redshift, confidence tests and comparisons with models 
through the remainder of Section 4. Finally, we summarize our results in Section 5. 

2. Ecirlier Studies 

Many past observations qualitatively support the statistical lensing picture. Character- 
istics of their respective studies, including galaxy and QSO catalogues, brightness ranges, 
QSO redshift ranges, as well as the observed association, can be found in Table 1. The 
'Detected' column indicates the type of statistically significant cross-correlations that were 
detected. '-I-', '0' and '— ' refer to positive-, no-, and ne^'atwe-correlations, respectively. The 
studies are divided into three groups from the top of the table to bottom, respectively: Top: 
Measurements of galaxy overdensities using incomplete QSO samples. Middle: Galaxy- 
QSO angular cross-correlations using complete sets of optically-selected (or in one case, UV- 
selected) QSOs. Bottom: Cross-correlations involving complete catalogues of radio-selected 
QSOs. Note: Romani & Maoz (1992) and Maoz (1995) are also mentioned in the introduc- 
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tory text, though they did not specifically measure overdensities or cross-correlations using 
a single range of QSO magnitudes and are as a result not listed on this table. 

Several early studies, including those found in the top portion of Table 1: Tyson (1986); 
Fugmann (1988, 1990); Hammer & Le Fevre (1990); Hintzen et al. (1991); Bartelmann 
& Schneider (1994), and Drinkwater et al. (1992); were conducted using heterogeneous 
or incomplete QSO samples, often involving portions of the Hewitt & Burbidge (1980) or 

Veron-Cetty & Veron (1989) QSO catalogues, and searched for associations of foreground 
galaxies within some angular scale around background QSOs relative to the field average. 
Not all of these studies necessarily attributed measured overdensities with gravitational 
lensing. However, even those studies that did not attribute them to lensing do, in retrospect, 
qualitatively agree with the magnification bias picture. 

While using optically-selected QSO catalogues, at least four studies found positive cor- 
relations, as expected on the basis of their specific QSO magnitude cuts. Webster et al. 
(1988) detected positive correlations between QSOs and galaxies on arcsecond scales and in- 
dicated that the foreground lensing mass must be much greater than that associated with the 
visible galaxies themselves. Rodrigues- Williams & Hogan (1994) observed an overdensity of 
LBQS QSOs with 1.4 < ^ < 2.2 and B < 18™5 within 6.6h-^ Mpc of Zwicky clusters (which 
have a typical redshift around 0.2) with 4.7(7 significance. Another study, by Williams & 
Irwin (1998) explored statistical correlations between LBQS QSOs and APM galaxies with 
a median redshift of 2; ^ 0.25. Significant positive correlations were found on angular scales 
< 100' , using a sample of QSOs at redshifts z > 1 with magnitudes Bj < 18™0. Finally, 
Gaztanaga (2003) used Sloan Digital Sky Survey Early Data Release galaxies and QSOs 
(albeit a heterogenous sample) and found cross-correlations on scales < 10'. Each of these 
studies measured cross-correlations that were a few times larger than what would be expected 
from CDM large scale structure lensing models. 

Even though the amplitudes of correlations are typically much greater than expected 
from models, the signal is most probably due to lensing because of the unique signature of 
lensing: Positive, negative and null correlations have been detected using QSO samples with 
the expected corresponding brightness cuts. This may hold true for some studies that were 
not originally looking into gravitational lensing as well. Boyle, Fong Sz Shanks (1988), for 
example, detected an 30% deficiency in the number oi B < 20^0 high-redshift UVX-selected 
QSOs at separations within 4' of foreground clusters of galaxies. Citing extinction due to 
dust as the simplest explanation, they ascribed an average extinction of Ay ^ 0"^2 to their 
set of clusters. However, an anti-correlation is expected within the magnification bias picture 
for QSOs with this limiting magnitude due to the resulting shallow slope, am- In fact. Groom 
& Shanks (1999) investigated the Boyle, Fong & Shanks (1988) anti-correlation under the 
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assumption that statistical lensing was the cause and found no discrepancy in interpretation 
between it and the positive correlation measured by Williams & Irwin (1998). Similarly, 
Romani & Maoz (1992) measured a ~ 25% dearth of high-redshift QSOs within 60' of nearby 
Abell Clusters and explained that it could be accounted for by the presence of ^4^ = 0™15 of 
extinction within the clusters. Later analysis by Maoz (1995) tested the presence of dust by 
looking for differences between the colors of QSOs that lay within 1° of Abell clusters with 
those outside of that separation. No difference was found and tight constraints on reddening 
were given: E{B — V) < 0^06 to 90% confidence. Just as in the case of Boyle, Fong & Shanks 
(1988), the anti-correlations that were measured in the Romani & Maoz (1992) work can 
hkely be explained in terms of magnification bias. Romani & Maoz (1992) show that the 
average QSO magnitude in their sample increases with decreasing QSO-cluster separation. In 
fact, the average magnitude reaches V ^ 19 for separations within 10'. Because the limiting 
magnitude is important in magnification bias, while neither individual QSO magnitudes 
nor limiting magnitudes of the QSO sample are given, it is quite likely that the slope, 
decreases with decreasing separation (and increasing limiting magnitude, inferred from the 
increase in the average magnitude). This would render a result consistent with weak lensing. 

A natural way to test for the influence of dust (Galactic or extragalactic) would be to 
compare cross-correlations involving optically-selected QSOs with those that utilize radio- 
selected ones. Many of the older studies hsted above do use radio-selected QSOs. In addition, 
several studies have incorporated radio-selected QSOs in order to determine whether extinc- 
tion is important. Comparing radio- and optically-selected subsets of QSOs from the Parkes 
and LBQS catalogues, respectively, with similar redshift distributions {z > 0.3), Benitez & 
Martinez-Gonzalez (1997) noted a ~ 99% significant excess of radio QSOs within 10' of 
COSMOS/UKST galaxies. The Parks excess is anti-correlated with the LBQS {Bj < 20°^5) 
sample at high confidence. And while magnification bias can account for part of the differ- 
ence, they argiie that dust could account for the rest if the extinction were present in clusters 
at a level consistent with Maoz (1995). Norman & Williams (2000) measured correlations 
between 1 Jy QSOs and APM galaxies on degree scales and found positive correlations that 
are inconsistent with a Galactic dust scenario while Benitez et al. (2001) used two complete 
samples of QSOs and found a similar result with the added caveat that the introduction 
of strong lensing could help to describe the discrepancy between lensing observations and 
models on smaller angular scales. 

In general, prior work in galaxy-QSO associations are at least qualitatively consistent 
with expectations from magnification bias, though some measurements yielded results that 
were significantly greater than predicted by CDM large scale structure lensing models. How- 
ever, inspection of Table 1 demonstrates that those studies using faint optical QSOs tend 
to have negative correlations while those that used bright QSOs found positive correlations. 
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This is entirely consistent with Equation 1. 

3. Data 

3.1. APM Catalogue 

Our galaxy data originated from the UKST SES Bj and R sky surveys, the photometry 
of which have been digitized using the Automated Plate Machine (APM) at the University 
of Cambridge (Irwin 1996). The APM galaxy catalogue, which has limiting magnitudes 
of Bj ^ 22.5 and R ^ 21, is complete to magnitudes of Bj ^ 20.5 and R ^ 19.5 with 
some plate-to-plate variation due to uncertainties in absolute photometric calibration and 
variations in sky background for each field. We have not sought to correct for the variation in 
photometric calibration, rather we adopt an analysis strategy that should remove its effect, 
as we will describe in Section 4.1. We culled objects that were detected and classified as 
galaxies in the R band with magnitudes in the range 18.5 < it! < 20 from 209 5?8 x 5?8 
UKST sky survey fields lying at declinations below -|-2.5°. 

Because the catalogue is derived from a photographic Schmidt survey, vignetting is 
a factor on individual UKST plate fields. It reduces the plate density of objects, but the 
magnitude of the effect is negligible within a radius of 2?7 from respective UKST field centers. 
Hence in order to avoid the contamination of the galaxy-QSO cross-correlation, we remove all 
galaxies lying more than 2? 7 from their respective field centers from our study. The locations 
of the galaxies used in our study, along with QSO positions can be found in Figure 1. 

We estimated the rcdshift distribution of galaxies in our magnitude cut using the fit given 
by Baugh & Efstathiou (1993). The parameters of this fit match the observed magnitude- 
binned galaxy redshift distributions from the Stromlo/APM redshift survey (Loveday et al. 
1992a,b), as well as the Durham/ Anglo- Australian Telescope faint galaxy redshift survey 
conducted by Broadhurst, Ellis & Shanks (1988) and the LDSS faint galaxy survey by CoUess 
et al. (1990, 1993). (Please see Section 4.1 of Maddox, Efstathiou & Sutherland (1996) for 
more details regarding the development of this fit.) The resulting distribution (with arbitrary 
normalization), shown as a sohd curve in Figure 2, peaks at a redshift of 2; 0.25 and falls 
off almost entirely by 2;~ 0.5 — 0.6. The galaxy sample does overlap the low-redshift portion 
of our QSO sample (sec Section 3.2), which we retain to investigate differences between the 
cross-correlation signals due to physical associations and magnification bias. 
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3.2. Hamburg-ESO Survey for Bright QSOs 

The Hamburg-ESO Survey for Bright QSOs is a survey that employs sUtless spectroscopy 
on the ESO Im Schmidt Telescope to find bright QSOs over an area that will eventually 
cover roughly 5000 0° of the southern sky. See Wisotzki et al. (1996) and Reimers, Koehler 
& Wisotzki (1996) for details regarding the methods and techniques used in the development 
of this survey. Currently, a preliminary portion (roughly half) of the final dataset is available 
to the pubfic. This subset contains 415 bright QSOs and Seyfert I nuclei with magnitudes 
in the range 13 < Sj < 17.5 and redshifts z < 3.2 that span an effective area of roughly 
3700 \J° (Wisotzki et al. 2000). Crowding of sources, especially on fields of low Galactic 
latitude, does result in the loss of spectra for roughly ~ 30% of the spectra on the fields. 
Nevertheless, the Hamburg-ESO survey QSO sample is particularly useful for this work 
because it is 99% spectroscopically complete over the well-defined ranges of magnitudes and 
redshifts given above and because it is devoid of redshift- and magnitude-related selection 
biases. The redshift distribution of Hamburg-ESO QSOs is shown as a dashed line histogram 
in Figure 2. 

Though Hamburg-ESO QSOs were identified on plates containing slitless spectra corre- 
sponding to ESO/SRC field position (not the UKST), QSOs selected for this study were those 
whose positions he within the 2?7 non- vignetted regions of the nearest respective UKST fields 
from which we had garnered our galaxy data. There is more overlap between UKST fields at 
high declinations than there is near the Celestial Equator. Therefore some high-declination 
QSOs lie within 2?7 degrees of multiple field centers. In these cases, we associated the QSO 
with the plate corresponding to the nearest field center. Those QSOs that lay further than 
2? 7 from any field center were not included in this study. 

The QSOs that fall within the apparent magnitude range found in the Hamburg-ESO 
Catalogue are typically intrinsically brighter than M* at most (especially at higher) red- 
shifts. This implies that the catalogue samples the steep portion of the QSO number- 
magnitude count distribution and that one would expect to measure positive galaxy-QSO 
cross-correlations d la Equation 1. 

4. Analysis 

4.1. Estimation of the Correlation Function 

The estimation of the galaxy-QSO cross-correlation is in many regards very similar to 
the determination of the galaxy auto-correlation, which has been conducted by many studies 
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in the past (Peebles 1980; Hewett, P. C. 1982; Maddox, Efstathiou & Sutherland 1996) 
and refinements have been incorporated in the technique (Hamilton 1993; Landy & Szalay 
1993), but there are some notable differences as well. 

Unlike autocorrelation estimates, the galaxy-QSO cross-correlation involves two sepa- 
rate catalogs, each with its own selection and window functions. In many cases, the selection 
function could also be an issue because one does not know a priori cither the underlying 
galaxy or quasar distribution from which survey samples are obtained. Furthermore, the 
window function of a survey could also influence the correlation estimates by introducing 
geometrical effects that conspire with the sky densities of objects to skew correlation esti- 
mates relative to the values one would obtain from fair samples that very well characterize 
the true distributions of galaxies and QSOs. 

An additional complication is the fact that the number of galaxies in our magnitude cuts 
from the APM catalogue is on the order of ~ 10*^ times greater than the number of QSOs 
in the Hamburg-ESO QSO catalogue. As a result, the uncertainty due to Poisson noise is 
greater relative to the signal for galaxy-QSO cross-correlations than it is for galaxy-galaxy 
auto-correlations . 

A number of autocorrelation estimators can be found in the literature, including: 

DD ■ (RR) ^ 

DD^2{DR) + {RR) 
^GGsie) jj^^ , (4) 

and 

Here, D and R refer to whether the galaxies were drawn from direct or randomly generated 
samples, respectively, in order to generate the pair counts DD, DR and RR. Angled brackets 
represent averages over a large number of realizations (100 realizations in the case of this 
work) of the total pair counts between combinations of random (artificial) galaxies. Within 
each realization, we use a number of randomly generated objects as there are real objects in 
order to preserve normalization (this is true in the case of cross-correlations as well). 

The most common estimator, uggi, has largely replaced the original one, ugga, in the 
literature and is often used in both auto- and cross-correlation (in a variant form) studies. 
This is because cugg^ is most severely affected by window effects of the correlators men- 
tioned here. It is a fact that has been known for quite some time, and has resulted in the 
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development of the newer correlation estimators, u>gg2 and ojggzi that Hamilton (1993) and 
Landy & Szalay (1993) developed in order to minimize the influence of contamination by 
the window and selection functions. These latter functions have the useful property that 
they are dependent on the galaxy-catalogue correlation function to higher orders than uiggi 
and should be less affected by selection and window effects as a result. 

We have adapted these autocorrelation estimators for use in our cross-correlation study: 

and, 

where the subscripts G and Q indicate whether the direct or random samples represent 
galaxies or QSOs, respectively. 

Like the galaxy autocorrelation estimators, each of the cross-correlation estimators are 
dependent to varying degrees on catalogue window functions and selection effects. Therefore, 
we have extended the formalism describing these dependencies found in Hamilton (1993) 
to the 2D auto- and cross-correlation estimators. Following the Hamilton (1993) formalism 
for a 2D catalogue, let us first define the catalogue selection function, $3., which could be a 
function of object brightness, color, position, etc., for catalogue, x, which would be either 
X = G 01 X = Q for galaxies or QSOs, respectively. Unlike Hamilton (1993), we will 
assume that the positional dependence of $ is 2-dimensional. If p is the surface density of 
objects within a catalogue, then the total number of objects within the catalogue, x, can be 
determined with 

Nobs,x = ^x^x, (10) 

so that if a catalogue contained all of the objects that lay within the bounds of the survey 
region, the selection function would be $ = 1. $ is technically a discrete function, but it 
is entirely possible to treat it as a continuous function, especially with large numbers of 
catalogue objects. 

If the catalogue window function is Wx-, then is defined so that 

= w^^x, (11) 
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and the cross-pair window is 

Wgq = wgq^g^q. (12) 

The pair weighting, wqq, need not be related to the point weightings, wq and wq. Then the 
direct pair counts will involve Nobs,x, weighted by a function, Wx, and the direct-direct pair 
counts for cross-correlation analysis will be given by 

DgDq = WGQNobs,GNobs,Q = WcQ^G^QUciyQ, (13) 

while direct-random pairs will be given by either 

DgRq = WGQNobs,G^Q = WgQ^g'^QI^G (14) 

or by 

RgDq = WGQNobs,Q^G = WGQ^Q^Gl^Q, (15) 

and finally the random-random pairs are given by 

RgRq = Wgq = wgq^g^q- (16) 

If Vg a-iid vq are the mean sky surface densities of galaxies and QSOs lying within 
respective catalogue windows, then the degree to which the represent the true full-sky average 
surface density, can be estimated with the sample overdensities gg and ctq, where 

Ox = — z— (17) 

for catalogue x. Fair samples that very well represent the true sky density of objects will 
have sample overdensities Ox ~ 0. Similarly, the object-catalogue correlation, ^x-, which is 
defined by either 

or 

should also be zero for a perfectly fair sample because the true object surface density in a 
sample, Vx-, would be equal to the true cosmic mean, v. 

Considering that the measured cross-correlation function will be. in general, a biased 
estimate of the true cross-correlation, we can define the windowed correlation. 
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which is equivalent to an unbiased cross-correlation in ideal circumstances. 

The results, giving the correlators and their influence by the catalogue sample overden- 
sities, aa,aQ, the catalogue- window correlations, (f)G,(pQ, and the object-catalogue correla- 
tions, logq, are shown in Table 2. Additional estimators (015 . . .wy) are also included. While 
these three estimators are not useful for determining the galaxy-QSO cross-correlations, they 
are helpful in gauging the degree to which randomly generated object distribiitions match 
real object distributions for modeling purposes. The most problematic terms in the third 
column of Table 2 are the first-order terms involving the catalogue-window correlations (the 
(p^s). It is then apparent that just like ujgga, the cross-correlation estimator, UJ4, is inferior 
because of its greater dependence on the catalogue-window correlations, whereas and 
are only dependent on them to second order. 

Direct pair counts were determined by counting all of the galaxies within angular sepa- 
ration bins 6^ ± Y around QSOs, using galaxies that were taken from one field per QSO only. 
In addition to eliminating the effect of vignetting, this eliminates the need to account for 
plate to plate photometric sensitivity variations across the APM catalogue. Furthermore, all 
pairs were drawn from the inner 2?7 portion of each UKST field. Though it does limit the 
angular scale over which we can estimate the cross-correlation to scales 9 < 5?4, this scale is 
rather large when compared to the scales over which prior studies have investigated and it 
prevents the contamination of the correlation signal by the plate to plate variations in the 
APM Catalogue photometric calibration and it serves to mitigate the infiuence of vignetting 
on the plate-density of galaxies. 

The pair counts involving randomly generated QSOs and galaxies were also treated in 
a manner that limited the inclusion of geometric effects. For each QSO in our sample, 100 
Circularly- Distributed Random (CDR) QSOs were generated that lay in a circle about the 
field center with a radius equal to the separation between the real QSO and the field center, 
and pair counts were averaged over the 100 realizations. Using CDR QSOs ensures that 
radial sensitivity gradients that may be present on the APM Schmidt plates cancel out, as 
will much of the edge effect introduced by the 2?7 radius boundary. Random galaxies for 
use in determining the (RgRq) and (RgDq) terms that are found in the estimators 002,^^3 
and UJ4 were created by randomly scattering points on the field so that the probability of a 
random galaxy being generated at some spot on the field as equal for all points within 2? 7 
of the field center and outside that radius. The number of randomly generated galaxies 
and the number of real galaxies on the corresponding UKST field were the same in every 
realization in order to preserve normalization. 
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4.2. Angular Cross-Correlations 

With combinations of pair counts in hand, we estimated the angular cross-correlations 
using each of the four estimators from the previous section with subsets of QSOs from a 
number of redshift bins of width Az = 0.2 in the range < z < 2. This is the first time 
that all four of these cross-correlation estimators have been compared using the same QSO 
and galaxy data sets and each estimator gives very similar results in each subset. A typical 
comparison of each of the estimators and the magnitude of their respective RMS uncertainties 
is shown in Figure 3 using QSOs within a redshift range of 0.8 - 1. As expected, the greatest 
discrepancies between them are at the smallest and largest galaxy-QSO separations where the 
variance (also shown) is highest. At small separations, the variance is due largely to Poisson 
fluctuations in the small number of galaxies located very near QSOs while the rise in variance 
at large pair separations is due both to increasing relative Poisson uncertainties from the 
increasing importance of edge effects near the limiting separation scale of 5?4 = 324' of our 
study use the dispersion among the individual QSO estimates as a conservative measure of 
uncertainty rather than the error in the mean due to the possible presence of the various 
systematic effects that were mentioned in the previous section. 

Among the individual correlators, ^4(6*), deviates the most from the others, especially 
at small and large separation scales where the stronger dependence of co^ on the window 
function becomes important. This is not surprising as UJ4 has long been considered an 
inferior estimator due to its relatively strong dependence on selection effects and window 
functions (Hamilton 1993). On the other hand, uji, U2, and cus all yield very similar results 
across, which means that Ui is apparently still an acceptable choice. To be sure, Landy 
& Szalay (1993) and Hamilton (1993) have shown that UJ2 and are less dependent on 
window and selection effects than ui, but they also require significantly more computational 
time in order to sample the additional random galaxy population that coi does not require. 
Hence if the cost of computation time exceeds the benefit of marginal signal improvement, 
then perhaps coi would prove useful. 

Furthermore, that each of the estimators, with their different dependences on the 
catalogue-object cross-correlations and relative over- or underdensity relative to the true 
object populations, are so similar imphes that to a good degree both the APM galaxy cata- 
logue and the Hamburg-ESO QSO catalogue are fair samples of their respective true object 
populations. This is significant because a number of earher studies (see Table 1) that have 
used the APM galaxy catalogue have found significantly stronger cross-correlations than ex- 
pected from models. If in some way, the APM did not represent a fair sample, statistical 
associations could conceivably result. 

Finally, the angular cross-correlation measurements are quite significant. While it is 
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true that in the case of a single QSO on an APM field, the galaxy counts in adjacent annular 
separation bins will be correlated at some level because of the clustering of galaxies on the 
sky and because of sensitivity variations on the plate, this effect is likely small because we 
have seen that each of the different estimators yield similar results despite varying sensitivity 
to the underlying galaxy distribution. And while individual separation bins are only detected 
to 2(7 on Figure 3, given that there are 415 QSOs in the Hamburg-ESO sample, which covers 
3700 n°, the average separation between QSOs is ~ 3°, so each of the small-separation angle 
bins should be nearly statistically independent because there should be little overlap of small- 
separation annuli around QSOs. As a result, considering that the small separation bins all 
show positive correlations and that they also follow an overall trend inconsistent with zero 
slope indicate that the detection of positive galaxy-QSO cross correlations is rather strong. 



If the cross-correlations discussed in Section 4.2 are due to magnification bias, then their 
amplitudes should monotonically rise with increasing redshift beyond the lensing plane (at 
z ~ 0.25) (Dolag & Bartelmann 1997; Bartelmann & Schneider 2001; Menard & Bartelmann 
2002). We integrated the angular cross-correlations over scales up to 30 and 60' in order 
to increase the signal to noise ratio enough so that we could divide our QSO sample into a 
number of subsets based on QSO redshift and explore the ^—dependence of the associations. 

Figure 4 shows the integrated cross-correlation vs. redshift. Because of low QSO num- 
bers per bin, we have used very wide (A^ = 0.4) redshift bins and we have therefore over- 
sampled to show the degree of variation with redshift. We do detect a positive signal on the 
order of 1 — 3% over all redshifts beyond the lens redshift, which is qualitatively consistent 
with magnification bias, and as expected, the Poisson noise is greater for the 30' plots than 
for the Q < 60' plot. This can be comphcated by the integral constraint. 



'"7 JQ 

Positive correlations at small separations must be balanced by negative correlations at higher 
separations. 

Since the correlators were found to be very similar in the previous section, one would 
expect the same of the integrated estimators in Figure 4. This is generally the case, although 
the relative differences between the integrated estimators is a bit greater than the relative 
differences between the estimators themselves. Because 



4.3. Cross- correlations vs. Source QSO Redshift 
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most of the signal contribution to a;(< 9o) comes from angular scales very near 9o- Perhaps 
the differences between the integrated correlation function is then due to the modification of 
the angular correlation profile by the integral constraint. Minor differences between angular 
correlations on scales ~ 9o would be exacerbated by the geometric scaling ~ 2tt9o- 

The overlap between our galaxy and QSO samples at low redshift (Figure 2) means 
that the low redshift {z < 0.4) signal in Figure 4 is due to physical associations rather than 
magnification bias. The 9 < 30' plot does in fact show a drop in the cross-correlation signal 
over redshifts in this range, along with a rise at redshifts greater than 0.4. This seems to 
indicate a fall off in the signal due to physical associations between galaxies and QSOs in 
front of the lensing plane, followed by an increase in the magnification bias signal involving 
QSOs beyond the lensing plane. The 9 < 60' plot does not show a significant dip between the 
physical and lensing associations. Perhaps this is due to a combination of the oversampling 
in redshift and due to the integral constraint. Separations of 60' are of the same order of 
magnitude as the 5?4 limiting scale of our study and the integral constraint may be altering 
the correlations on these scales. 



4.4. Confidence Tests 

Individual angular cross-correlation functions from redshift bins with width Az — 0.4 
tend to be a bit noisy. Therefore, it is imperative to create a significance test for the 
integrated correlation function. Working toward this end, we have developed a Monte Carlo 
algorithm that generates a large number of realizations of the integrated cross-correlation 
measurements in the same manner as those from the last section. It utilizes our APM 
galaxy sample along with a large number of artificially generated QSO catalogues in order 
to determine the significance of the integrated correlation funtion measurements from the 
last section. 

The algorithm works as follows: First, for each QSO in the Hamburg-ESO Catalogue, 
we generated one random QSO so that the probability lying at a particular point within 
2?7 of the field center for this first type of QSO was equal for all points and outside of 
that radius (as we did for our randomly generated galaxies, see Section 4.1). We shall call 
these Uniformily- Distributed Random (UDR) QSOs. Each UDR point was located on the 
same field as the real QSO that it replaced. Furthermore, in order to preserve the number 
of QSOs per redshift and magnitude bin, each was also attributed the same magnitude and 
redshift as the replaced real QSO. Next, for each UDR QSO, we generated 100 CDR QSOs 
(sec Section 4.1) whose positions were distributed along a circle about the field center with a 
radius equal to the distance between the UDR QSO and the field center, as in our real-QSO 



- 15 - 



analysis. Finally, we determined the correlation function,a;i, from the resulting pair counts. 
This process was repeated 1000 times, effectively creating 1000 random realizations of the out 
cross-correlations. We interpret the distribution in values of the artificial cross-correlation 
functions to be indicative of the magnitude of the uncertainties in our measurement of the 
cross-correlation function. 

Figure 4 shows the 50, 75, and 99% confidence limits that we derive from this interpre- 
tation, corresponding to the ranking of the values of 1000 cross-correlation realizations in 
each redshift bin. As the diagram indicates, the physical association signal is detected to 
high significance, especially in the lowest redshift bin, corresponding to source QSOs with 
redshifts < 0.4. The magnification bias signal seems to be most significantly detected in the 
9 < 30' frame at redshifts around 0.8 and 1.3, although the overall detection is somewhat 
marginal - with confidence levels on the order of 75%. 

We do consider this Monte Carlo confidence analysis to be rather conservative, however. 
This is due to the fact that randomly generated QSOs can be either positively or negatively 
associated with APM galaxies and the confidence contours illustrate the ranking of individ- 
ual realizations. The 99% contour therefore samples the very extremes of the underlying 
distribution, which seems to have a rather wide tail. Hence we believe the confidence limits 
are somewhat underestimated. 

It is also useful to compare the galaxy-QSO results with those of cross-correlations 
between objects for which there should be no net cross-correlation. Therefore, we selected 
objects hsted as stars in the APM catalogue with magnitudes 16 < i? < 17.5. Aside from 
the contamination of our stellar sample by galaxies on the faint side of our magnitude and 
object class cut, there is no reason that Galactic stars should be correlated in any way with 
QSOs that lie at cosmological distances. Figure 5 shows the strength of the integrated Star- 
QSO cross-correlation vs. redshift that was obtained in the same manner as the galaxy-QSO 
cross-correlation. No real dependence on redshift can be seen as each of the cross- correlation 
estimators oscillate above and below in a manner that would be expected for random noise. 
Indeed, the star-QSO cross-correlation can serve to provide an independent indication of the 
noise level for the galaxy-QSO signal that is consistent with our Monte Carlo confidence 
estimate. 



4.5. Modeling 

We modeled the signal in Figure 4 based on a hnear CDM large scale structure lensing 
model, given by Bartelmann & Schneider (2001), which gives the expected correlation 
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function from an assumed mass power spectrum. This particular model is a simplified version 
of that given by Dolag & Bartelmann (1997) in which an approximation is made by assuming 
that all galaxies are found at their median redshift. This approximation incurs errors on the 
order of 10%. The mass power spectrum used for our model was Case 1 obtained from the 
Appendix G of Bardeen et al. (1986) (BBKS): A CDM power spectrum with adiabatic 
fluctuations. We assumed a value of the present-day mass fluctuation normalization of 
CTg = 1, a matter density parameter of ftm — 0.3, cosmological constant Q\ — 0.7, and 
expansion rate h — 0.75, resulting in a shape parameter of our modeled mass power spectrum 
was r = 0.225. The power spectrum at the lens plane redshift was then determined by scahng 
the present-day power spectrum with a factor of (1 -|- ziens)~^- 

Using our mass power spectrum and the assumption that all of the galaxies in our sample 
lie on a plane with a redshift z — 0.25, we determined the expected strength of the angular 
cross-correlations using the Bartelmann & Schneider (2001) approximation while varying 
the source redshift over the range represented by our QSO sample. Each of the correlation 
signals were integrated over separations up to 30 and 60', respectively. The model predicts an 
increase in the strength of the signal with source distance, provided that the sources lie behind 
the lensing plane. This can be seen in Figure 4, which agrees well qualitatively with the 
measured signal from the combination of Hamburg-ESO QSOs and APM galaxies. However, 
consistent with a large number of other studies (see the Introduction) , the integrated modeled 
correlations were significantly lower (roughly a factor of 20) than our measured correlations. 
We will explore this further in the Discussion section. 



4.6. Correlation scales 

As wc have mentioned earlier, we see in Figure 4 that QSOs with redshifts z < 0.4 are 
physically associated with the APM galaxy sample, while those at higher redshift are likely 
affected by magnification bias. We have separated our QSO sample into two redshift ranges 
in Figure 6, z < 0.2 and 0.6<2;<1.0, in order to see whether the differing assumed physical 
mechanisms result in different correlation functions. Outside of differences on angular scales 
greater than 130', the low and high redshift sample have very similar profiles. Both possess 
'correlation lengths' (defined here as 9o such that uj{9o) = 0, as opposed to the standard 
definition of the physical length, Tq, where ^(ro) = 1) of ~ 110', which corresponds to scales 
of roughly 23 Mpc at the redshift of our APM galaxies {z ~ 0.25) (assuming h = 0.75). 
These scales are much larger than the typical sizes of groups and clusters and they seem to 
indicate that both physical and lensing associations are due to Large Scale Structure. That 
the physical and lensing associations have roughly the same correlation length may be due 
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to the fact that the same structure is hkely responsible for both signals. 

4.7. Galactic Dust 

With such a large disagreement between the positive cross-correlations that are observed 
in this study with model expectations, it is appropriate to investigate whether these cross- 
correlations could in fact be generated by Galactic dust. If the QSOs in the Hamburg-ESO 
sample happen to lie in portions of the sky that are relatively unobscured by intervening 

Galactic dust, it could conceivably explain the amplitude discrepancy. Rodrigucs- Williams & 
Hogan (1994) showed that the overdensity due to differences in levels of extinction, SAv{0), 
between two lines of sight separated by an angle, 9, is given by 

q{e) = io-2-5-"-5^W ~ 1 _ 5.76 • amSA{e). (23) 

If dust were the sole cause of the observed cross-correlations, then the maximal overdensity 
produced by dust would be roughly proportional to the differential extinction. 

The Rodrigues- Williams & Hogan (1994) analysis determined the extremes in extinction 
across LBQS fields using the Burstein & Heiles (1982) maps and found that the observed 
galaxy-QSO cross-correlations between LBQS QSOs and APM galaxies were in fact anti- 
correlated with the cross-correlations that would have been expected if they were due solely 
to Galactic extinction. 

We used a slightly different technique in which we interrogated the Schlcgcl, Finkbeiner 
& Davis (1998) Galactic E{B — V) reddening maps using the IDL dust_getval routine, 
which is provided by Schlegel & Finkbeiner (1999), to derive and compare the differences 
in extinction levels between Hamburg-ESO QSO positions and large-scale averages within 
surrounding annuli with inner and outer radii of 135' and 150', respectively. We define the 
differential extinction by 

SAvie) = Av{< 15') - Av{lS5 - 150'). (24) 

The Schlegel, Finkbeiner & Davis (1998) reddening map has a resolution of 6™1 x 6™1 
and local reddening values were interpolated from the adjacent surrounding pixels. We 
determined extinction levels from the Schlegel, Finkbeiner & Davis (1998) reddening map 
by adopting a single value for the ratio of total to selective extinction Ry — 3.1 for every 
Hamburg-ESO QSO. 

Meanwhile, we determined the galaxy number densities from the same regions and 
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measured the relative overdensities between the QSO positions and the field: 




A^(135 - 150') ■ 



(25) 



The reason for the specific choice of the annular region size was that it represents the ap- 
proximate angular scale at which the galaxy-QSO cross-correlations, q{6) — 1 = uj{9) are 
near zero for our sample. Figure 7 shows that the relative overdensities and the differential 
extinctions between the QSO positions and the field are uncorrelated. In addition, the solid, 
dotted and dashed lines in Figure 7 show the maximal overdensities that would be expected 
from dust extinction alone, given slopes of = 0.6,0.4, and 1.0, respectively. The slope of 
0.6 corresponds roughly to the slope of the QSOs in our sample, while 0.4 corresponds to the 
critical slope at which there would be no expected lensing signal. If dust were the sole cause 
of the cross-correlation signal, one would expect that the magnitude of the overdensities 
(underdensities) would lie somewhere between a value of ^ = 1 and the curve derived using 
the appropriate QSO sample slope. Yet, the most of the QSO overdensities fall outside this 
bound. 



An important result of this work is the detection of a net positive correlation signal 
over the entire redshift range for which we have significant numbers of QSOs (see Figures 3 
and 6). Though our Monte Carlo simulations suggest a modest level of confidence for the 
integrated correlations, our having detected a consistently positive signal over all redshift 
intervals suggests that perhaps a somewhat higher level of confidence could be placed on 
the result. The Monte Carlo results also belie the fact that the angular cross-correlations 
from various redshift intervals arc quite significant and that they indicate strong detections 
of positive galaxy-QSO cross-correlations. We anticipate that when the final Hamburg-ESO 
QSO Catalogue is finished, effectively doubling the size of the current sample, the total S/N 
ratio should in principle improve by a factor of roughly -\/2, which would lend a significantly 
higher level confidence in the cross-correlation measurements. 

The measurements of statistical associations between APM galaxies and Hamburg-ESO 
QSOs occur on angular scales 9 < 110'. The overall shape of the correlations for physical 
and lensed associations are very similar. Meanwhile, the correlation length (equivalent for 
both types of associations, corresponding to ~ 20Mpc scales at ziens ~ 0.25 suggests that 
large scale structure is doing the lensing rather than more compact objects such as galaxy 
clusters. 



5. 



Discussion and Conclusions 



Despite the large order of magnitude discrepancy between our modeled and observed 
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integrated cross-correlation functions, we are consistent with several other studies, including 
Beni'tez & Martmez-Gonzalez (1997); Williams & Irwin (1998); Norman & Williams (2000) 
and Benitez et al. (2001), because they have also noted such large discrepancies. A number 
of potential causes for the mismatches have been discussed in the literature. On small 
angular scales less than a few arcminutes, correlations can be explained by the halo model, 
which indicates that the strength of the signal is very dependent on the halo occupation 
properties of various types of galaxies. Depending on the specific types of galaxies used, the 
halo model predicts an enhancement of a factor of 2 - 10 on small angular scales relative 
to nonhnear clustering models by Bartelmann (1995) and Dolag & Bartelmann (1997) 
(Jain, Scranton & Sheth 2003). There is still a discrepancy on the > 30' scales within this 
study, however, and a number of other phenomena have been investigated to explain this. 
Williams (2000), for example, investigated extending lensing formalism by including higher 
order general relativistic terms, but this results only in a ~ 10% increase in the theoretical 
cross-correlation estimate. Many other studies have suggested dust, but extragalactic dust 
would tend to reduce cross-correlation, but there is probably little or no dust in clusters 
(NoUenberg et al. 2003). Galactic dust could create positive correlations, but one would 
also expect to see positive correlations between faint stars and QSOs, which we do not detect. 
Studies involving radio-selected QSOs (see Table 1), which should be unaffected by dust, are 
also inconsistent with significant amounts of Galactic obscuration. We have also compared 
various cross-correlation estimators. Each are dependent to varying degrees on selection 
effects and window functions, but each yield very similar results. This likely implies that at 
least the APM and Hamburg-ESO catalogues are fair samples of their respective populations 
and that systematic effects are not inflating the cross-correlation estimates. 

So in order to further investigate the discrepancy between the model and the obser- 
vations, we have also measured the slope of the QSO number- magnitude counts, as a 
function of redshift using the same sample of QSOs for each redshift bin that was described 
in Section 4.3 (the binwidth is, again, = 0.4). Figure 8 shows that is indeed above 
the critical value of 0.4 for all redshifts, with perhaps a shght overall increase in the slope 
with increasing QSO redshift. Typical values are on the order of am ~ 0.6. The error bars 
represent the la RMS variation in the slope of Az — 0.1 subsets of the Az — 0.4 bins of 
our QSO sample and uncertainties at redshifts z > 1.5 are underestimated due to low QSO 
numbers. 

Using the measured am with the measured correlation strengths, we determined the 
degree of magnification present, assuming that the signal is due to magnification bias (see 
Equation 1). The magnification is plotted vs. measured correlation strength in Figure 9, 
using the same data found in Figures 4 and 8. Because the physical associations do not obey 
Equation 1, the estimated magnification should not necessarily compare with that expected 
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from the magnification bias affect points. Therefore, we are able to separate those points 
corresponding to low-redshift physical associations from the lensed points at higher redshift 
with this plot. 

Finally, because of the rather large amplitude discrepancy with the models, we investi- 
gated whether the observed overdensities could be due to extinction by Galactic dust. As 

mentioned in Section 4.7, there is no correlation between the measured galaxy overdensities 
about the QSOs and the differences in extinction level between the QSO lines of sight and 
the field. If dust were the cause, one would have difficulty in explaining the observed un- 
derdensities that have been measured in other studies that employ faint QSOs, as expected 
from magnification bias. Patchy Galactic dust can only produce positive cross-correlations, 
not anti-correlations. Both of these hues of evidence suggest that Galactic dust is not the 
cause of the large amplitude in the cross-correlation signal. Furthermore considering that 
the signal appears to behave as expected, qualitatively, for magnification bias, it seems likely 
that this is the cause of the observed galaxy-QSO cross-correlations. 

Yet we are still left with the long-standing problem of the discrepancy with models. 
There remain a number of avenues to explore in hopes of solving this problem. Dust, whether 
Galactic or extragalactic in origin, are unlikely to skew the values of the cross-correlations 
significantly. Galactic dust is incapable of explaining anti-correlations and extragalactic dust 
can not result in positive correlations. Both types of correlations are observed in accordance 
to magnification bias. We will investigate the effect of catalogue masks and selection effects 
in Nollenberg & Williams (2005), though there are still a number of other places to look for a 
solution. For example, perhaps models need to account for multiple or nonlinear deflections 
for a given light ray as it travels from source to observer. There has also been a suggestion 
by Bassett and Kunz that since galaxy-QSO cross-correlations are one of two types of studies 
(the other involving Type-la supernovae) in which both the luminosity and angular diameter 
distances are sampled, there may be an asymmetry between the two distances (Bassett & 
Kunz 2003, 2004). This asymmetry could in principle be caused by essentially any process 
that would hinder a photon from reaching an observer, including Compton scattering off the 
ionized IGM or something much more exotic such as photon decay. 

We thank Lutz Wisotzki for his provision of the initial release of the Hamburg-ESO 
QSO catalogue and for his very useful comments regarding the preparation of the survey 
and his equally helpful suggestions for this paper. 
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Table 2: 2D Estimator dependence on catalogue characteristics 



Correlator 
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Fig. 1. — Sky positions of APM galaxies (grey points) and Hamburg-ESO QSOs (black 
asterisks) used in this study. Left: Equatorial coordinates. Right: Galactic coordinates. 
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Fig. 2. — Redshift distributions of Galaxies and QSOs used in this study. The histogram 
represents the QSOs in the Hamburg-ESO QSO sample while the solid line shows the esti- 
mated redshift distribution (with arbitrary normalization) for our sample of roughly 4.5-10^ 
APM Galaxies with magnitudes 18.5 < -R < 20, following Maddox, Efstathiou & Sutherland 
(1996). 
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Fig. 3. — A comparison between correlation function estimators, using a sample of QSOs 
with redshifts O.S < z < 1.0 and magnitudes Bj < 17.5. Top: Angular correlation function 

vs. pair separation angle. The correlators a;2 and C1J3 are very similar in this case. Bottom: 
RMS Uncertainty in the correlation function calculated using the distribution of correlation 
functions of individual QSOs within the sample. The errorbars in the top panel show the 
RMS uncertainty for ui. 
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Fig. 4. — Measured and modeled integrated correlation functions vs. redshift with confi- 
dence limits. Each of the four correlation estimators from the Methods section are shown. 
Top: Correlation functions integrated over separations < 60'. Bottom: The same, but now 
integrated over separations < 30'. The bold dark line in each case represents a large scale 
structure lensing model with a linear CDM mass power spectrum, following Bartelmann & 
Schneider (2001). Dotted, dashed and long-dashed lines represent 50, 75 and 99% con- 
fidence limits from 1000 iterations of our monte carlo error estimation. Note the positive 
signal detected at redshifts z > 0.25 and that cus nearly coincides with coi on both panels. 
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Fig. 5. — Same as Figure 4 only substituting galaxies with a selection of APM stars with 
magnitudes between 16 < i? < 17.5. One would expect no correlation between stars and 
distant QSOs. Unlike the positive signal seen in Figure 4, we indeed find no correlation 
between stars and QSOs. Offsets between the various estimators are due to differing shapes 
of the individual correlation functions from each estimator. Some contamination of the 
stellar sample by galaxies is possible. Several estimators are in close agreement for the 60' 
correlations while 002 and uj^ are most similar in 30' correlations. 
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Fig. 6. — Correlation functions from two redshift cuts. Empty Squares: All 113 QSOs with 
< ^ < 0.2 (physical associations). Solidd dots: The 79 QSOs with 0.6 < 2; < 1.0, corre- 
sponding to the location of the strongest signal in the integrated correlations from Figure 4 
(magnification bias). Despite their different physical cause, both correlation functions appear 
to be very similar considering the statistical uncertainties. 
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Fig. 7. — The observed overdensities of galaxies within 15' of QSOs with respect to galax- 
ies between 135 — 150' is plotted here versus the difference in Galactic extinction between 
the same separation scales. The points are divided between high- (60 - 90°), mid- (30 - 
60°), and low-Galactic Latitude (0 - 30°) samples (empty circles, stars, and filled circles, re- 
spectively). The points appear completely uncorrelated. In addition the expected maximal 
overdensities are plotted as a function of differential extinction for cases in which the QSO 
cumulative number-magnitude count slope — 0.6 (solid), 0.4 (dotted), and 1.0 (dashed). 
(See Equation 23.) 
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Fig. 8. — Average slope of the QSO number-magnitude counts plotted vs.redshift. Redshift 
bins have width Az = 0.4 and each bin contains the same QSOs as the bins found in Figure 4. 
The value of the slope of the QSOs within the Hamburg-ESO Catalogue is > 0.4 for all 
redshifts, implying the expectation of positive correlations due to magnification bias. Error 
bars represent the la RMS of slopes from subsets of the data. Uncertainties given at redshifts 
between 1.5 < z < 2.0 are artificially low due to low QSO numbers per bin. 
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Fig. 9. — The required magnification to generate the measured integrated correla- 
tion, assuming lensing and given the measured slope of the QSO number-magnitude counts, 
am{z) for each redshift bin from Figure 4. X's, empty squares and filled squares correspond 
to points from uj{< 15'),a;(< 30') and a;(< 60'), respectively. The points can be separated 
between the two lowest redshift bins, and the rest. (The dotted line separating the two pop- 
ulations is, outside the fact that it divides the two populations, arbitrary.) The requirement 
to see such large magnifications to account for the measured correlations when compared to 
the rest of the points implies a different mechanism is the cause of the measured correlations: 
Physical associations. 



